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Abstract 


A search for new particles that decay into top quark pairs is reported. The search is per¬ 
formed with the ATLAS experiment at the LHC using an integrated luminosity of 20.3 fb“^ 
of proton-proton collision data collected at a centre-of-mass energy of = 8 TeV. The 
lepton-plus-jets final state is used, where the top pair decays to W^bW~b, with one W boson 
decaying leptonically and the other hadronically. The invariant mass spectrum of top quark 
pairs is examined for local excesses or deficits that are inconsistent with the Standard Model 
predictions. No evidence for a top quark pair resonance is found, and 95% confidence-level 
limits on the production rate are determined for massive states in benchmark models. The 
upper limits on the cross-section times branching ratio of a narrow Z' boson decaying to top 
pairs range from 4.2 pb to 0.03 pb for resonance masses from 0.4 TeV to 3.0 TeV. A narrow 
leptophobic topcolour Z' boson with mass below 1.8 TeV is excluded. Upper limits are set on 
the cross-section times branching ratio for a broad colour-octet resonance with T/m = 15% 
decaying to tt. These range from 4.8 pb to 0.09 pb for masses from 0.4TeV to 3.0TeV. A 
Kaluza-Klein excitation of the gluon in a Randall-Sundrum model is excluded for masses 
below 2.2 TeV. 
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1 Introduction 


Many models of physics beyond the Standard Model (SM) predict production at the LHC of additional 
particles with masses near the TeV scale. This paper presents a search for such heavy particles decaying 
to top quark pairs {tt) using data from proton-proton collisions collected at a centre-of-mass energy 
= 8 TeV with the ATLAS detector, corresponding to an integrated luminosity of 20.3 fb“^ 

Searches for production of heavy particles that decay to tt are of high interest at the LHC due to the 
role that the top quark plays in many models of physics beyond the SM (BSM). The top quark is the 
most massive of the fundamental particles in the SM, and it can have a large coupling to heavy Higgs 
bosons. Thus, heavy Higgs bosons in two-Higgs-doublet models [1,2] can have a large branching ratio to 
tt final states. Furthermore, many models that propose alternative mechanisms for electroweak symmetry 
breaking (EWSB) incorporate new heavy particles with a larger coupling to tt than to lighter quarks. 
Examples include strong EWSB models such as topcolour-assisted technicolour [3] (TC2) and Composite 
Higgs [4-11] scenarios. Models with warped extra dimensions [12-15] form an additional class of models 
that predict heavy particles that decay to W pairs. In such models, the heavy particles are the counterparts 
of the gluon and graviton. 

The search starts by selecting events with one isolated charged lepton (electron or muon), missing trans¬ 
verse momentum (whose magnitude is denoted by and hadronic jets, which are compatible with 

tt W'^bW~b, with one W boson decaying leptonically and the other hadronically. At least one of the 
hadronic jets is required to be consistent with having originated from a (7-quark. An estimator of the tt 
invariant mass (mjp”) is constructed with the events divided into two orthogonal classes by topology: 
the boosted topology, where the decay products of the hadronically decaying top quark are expected to 
be fully enclosed within one large-radius jet, and the resolved topology where four small-radius jets are 
reconstructed and attributed to the bbqq' quarks. The spectrum is scanned for localised deviations 
relative to the expectations from the background processes. The compatibility of the data and expectations 
is assessed, and limits on the production cross-section for new particles are set. 

This search is designed to be sensitive to the production of any new particle that decays to tt. Nonetheless 
the selection efficiency and acceptance can differ befween particular model choices. Hence fhe sensitivity 
fo a variety of differenf new particles was evaluafed fo quanlify fhe performance of fhe search. The 
benchmark models adopfed in fhis search include colour-single! and colour-ocfef bosons wifh spin 0, 1 
and 2 and masses from 0.4 fo 3 TeV. The resonance widfh for fhe specific models varies from very narrow 
(1%) fo a size similar fo thaf of fhe experimenfal resolution (15%). Eurfhermore, fhe dependence of fhe 
limifs on fhe resonance widfh is explored for heavy gluons up fo a widfh of 40%. Wifh fhese resulfs, if 
is possible fo inferpref fhe cross-secfion limifs in fhe confexf of a search for ofher new particles wifh fhe 
same producfion modes. 

In addition fo making use of 8 TeV dafa, fhis search incorporafes several improvemenfs wifh respecf to 
the previous most sensitive ATEAS collaboration search for the same signature [16], which used a similar 
strategy. Trimming [17] is now used to mitigate pile-up effects on large-radius jets, a new procedure is 
applied to recover efficiency loss from removal of overlap befween election and jef objecfs in fhe defecfor, 
and variable is exploifed to remove events from non-ti background processes. Eurthermore, the limits 
on production of spin 2 bosons and on heavy gluons of varying width represent an extension compared to 
previous searches by the ATEAS and CMS [18] Collaborations. 
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2 Models tested 


The details of the benchmark models considered in this search are reviewed below. Interference between 
these processes and SM tt production is not considered in this search. 


2.1 Spin-1 colour singlet 

The first class of models explored produces spin-1 colour-singlet vector bosons, Z'. This search uses 
topcolour-assisted technicolour [3, 19,20] as a benchmark. This is a leptophobic boson, with coup¬ 
lings only to first- and third-generation quarks, referred to as Model IV in ref. [19]. The properties of 
the boson are controlled by three parameters: cot On, which controls the width and the production cross- 
section, and /i and / 2 , which are related to the coupling to up-type and down-type quarks respectively. 
Here f\ = I and /2 = 0, which maximises the fraction of that decay to tt. The parameter cot On 
is tuned for each mass point such that the resonance has a width of 1.2% of its mass. To account for 
higher-order contributions to the cross-section, the leading-order calculation is multiplied by a factor of 
1.3 based on calculations performed at next-to-leading order (NLO) in QCD [21,22]. 

Constraints on have been set by the CDF [23,24] and DO [25] collaborations using data from proton- 
antiproton collisions at the Tevatron. Previous constraints at the LHC were set using proton-proton 
collisions at = 7 TeV with an integrated luminosity of 5 fb“^ by the ATLAS [16,26] and CMS [27-29] 
collaborations, and using 20fb“' of = 8 TeV data by the CMS Collaboration [18]. For narrow (wide) 
^TC 2 width 1.2% (10%) the strongest lower bound on the allowed mass is 2.1 TeV (2.7 TeV) from the 
search performed by CMS at = 8 TeV. 

2.2 Spin-1 colour octet 

The second class of models considered produces spin-1 colour-octet vector bosons. Specifically, heavy 
Kaluza-Klein gluons, q'kk. as produced in Randall-S undrum (RS) models with a single warped extra 
dimension [30], are used as a benchmark in this search. In this model, the ^kk has a nominal width of 
15.3% of its mass. Previous searches using ^[s = 1 TeV ATLAS data [16] exclude a ^kk with a mass 
less than 2.1 TeV. The CMS Collaboration searched for similar resonances [18], using a slightly different 
benchmark model [31]. The CMS choice leads to a larger natural width of 20% and a larger production 
cross-section, and, for such a scenario, CMS excludes the existence of a ^kk with mass less than 2.5 TeV. 
In the analysis presented here, the sensitivity to the width of the colour octet is also tested for widths from 
10% to 40% of the resonance mass. 


2.3 Spin-2 colour singlet 

The third class of models explored in this search produces spin-2 colour singlets, such as Kaluza-Klein 
excitations of the graviton, Gkk- The search uses a Randall-Sundrum model with extra dimensions where 
the SM fields are in the warped bulk and the fermions are localised appropriately to explain the flavour 
structure of the SM [12,32,33]. This kind of graviton is commonly referred to as a “Bulk” RS graviton 
and is characterised by a dimensionless coupling constant k/Mpi ~ 1, where k is the curvature of the 
warped extra dimension and Mpi = Mpi/ is the reduced Planck mass. For such gravitons, decays 
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to light fermions are suppressed, and the branching ratio to photons is negligible. The branching ratios 
to tt, WW, ZZ and HH are significant. In the model used, klM^x is chosen to be 1, and the Gkk width 
varies from 3% to 6% in the mass range 400-2000 GeV. The branching ratio of Gkk decay into a tt 
pair rapidly increases from 18% to 50% between 400 and 600 GeV, plateauing at 68% for masses larger 
than 1 TeV. There have been no previous direct searches for such gravitons in the tt decay channel. The 
ATLAS Collaboration used the same model to explore the Gkk —^ channel [34] and excluded Bulk 

RS Gkk with mass less than 740 GeV. The CMS Collaboration performed searches in the Gkk ^ 2Z and 
Gkk —^ WW decay channels [35] but did not consider the case of Bulk RS gravitons with klMpi > 0.5. 

2.4 Spin-0 colour singlet 

The last class of models examined here produces colour-singlet scalar particles via gluon fusion which 
decay to tt. The approach previously adopted by the CMS Collaboration [18] is followed, in which narrow 
scalar resonance benchmarks are generated while the interference with SM tt production is neglected. 
Even though such signals with negligible interference are not predicted by any particular BSM model, 
they can be used to evaluate the experimental sensitivities and set upper limits on the production cross- 
sections. The CMS Collaboration excluded such resonances with production cross-sections greater than 
0.8 pb and 0.3 pb for masses of 500 and 750 GeV, respectively. 


3 The ATLAS detector 

The ATLAS experiment [36] is a multipurpose particle physics detector with forward-backward sym¬ 
metric cylindrical geometry.^ The inner detector (ID) consists of multiple layers of silicon pixel and 
microstrip detectors and a straw-tube transition radiation tracker and covers a pseudorapidity range of 
I 77 I < 2.5. The ID is surrounded by a superconducting solenoid that provides a 2 T axial magnetic 
held. The calorimeter system, surrounding the ID and the solenoid, covers the pseudorapidity range 
\r]\ < 4.9. It consists of high-granularity lead and liquid-argon (LAr) electromagnetic calorimeters, a steel 
and scintillator-tile hadronic calorimeter covering \t]\ <1.7 and two copper/LAr hadronic endcap calor¬ 
imeters covering 1.5 < I 77 I < 3.2. Forward copper/LAr and tungsten/LAr calorimeter modules complete 
the solid-angle coverage out to \t]\ = 4.9. The muon spectrometer resides outside the calorimeters. It 
consists of multiple layers of trigger and tracking chambers within a system of air-core toroids, which 
enables an independent, precise measurement of muon track momenta for \r]\ < 2.7. The muon trigger 
chambers cover \t]\ < 2.4. 


4 Data and Monte Carlo samples 

This search is performed in proton-proton collision data at = 8 TeV collected with the ATLAS 
detector in 2012. The data are only used if they were recorded during stable beam conditions and with 


’ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the !/-axis points 
upward. Cylindrical coordinates (r, (p) are used in the transverse plane, (p being the azimuthal angle around the beam pipe. 
The pseudorapidity is defined in terms of the polar angle 9 d& rj = -lntan(6/2). The distance in r]-(p space is commonly 
referred to as AI? = + (A;;)^. 
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all relevant subdetector systems operational. Lepton-plus-jets events were collected using single-electron 
and single-muon triggers with thresholds chosen in each case such that the efficiency is uniform for 
leptons satisfying offline selections including transverse momentum pj > 25 GeV. The ATLAS muon 
trigger system suffers from a 20% inefficiency, relative to the offline event selection used in this analysis, 
largely due to a lack of geometrical coverage by muon chambers owing to support structures in those 
regions [37]. To mitigate this loss of efficiency, a large-radius-jet trigger, which triggers on anti-k( jets 
(see section 5) with radius parameter R = 1.0, was used to collect muon-plus-jets events which failed 
the muon trigger. This large-radius-jet trigger recorded 17.4 fb“^ data. The chosen trigger threshold 
yields a uniform efficiency exceeding 99% for events containing a large-radius jet with reconstructed 
Pt > 380 GeV. For tt events with invariant masses above 1.5 TeV, this addition increased the overall 
trigger efficiency in the muon channel to 96%^. 

Simulated Monte Carlo (MC) samples were used for signal processes, as well as background processes 
producing jets and prompt leptons. The MC samples are employed to develop the event selection, provide 
SM background estimates, and evaluate signal efficiencies. Background contributions from processes in 
which there are no genuine prompt isolated leptons were estimated directly from the data as described in 
section 7. The MC samples were processed through the full ATLAS detector simulation [38] based on 
Geant4 [39] or through a faster simulation making use of parameterised showers in the calorimeters [40]. 
Additional simulated proton-proton collisions generated using Pythia v8.1 [41] were overlaid to simulate 
the effects of additional collisions from the same and nearby bunch crossings (pile-up). All simulated 
events were then processed using the same reconstruction algorithms and analysis chain as the data. The 
simulated trigger and selection efficiencies were corrected to agree with the performance observed in 
data. 

Production of Bulk RS gluon and graviton signals was modelled using MadGraph5 [42] interfaced with 
Pythia v8.L For the gluon, the MSTW2008LO parton distribution function (PDF) set [43] was used, 
while for the graviton, the CTEQ6L1 [44] PDF set was used. The Z' signal was modelled using Py¬ 
thia v8.1 with the MSTW2008LO PDF set. Heavy scalar signal samples were generated using Mad- 
Graph_aMC@NLO [45] with LO matrix elements and the CTEQ6L1 PDE set. 

Pair production of top quarks is the dominant background in this search. It was simulated using the 
Powheg-Box [46-49] generator r2330.3 interfaced with Pythia v6.427 [50] with the Perugia 2011C [51] 
tune and the CTIO [52] next-to-leading-order PDE set. In the MC generation of Powheg-Box events, 
the parameter hdamp was set to the top quark mass, corresponding to a damping of high-px radiation, 
in order to achieve good agreement with the differential cross-section measurements [53]. Alternative 
samples of tt events, used to evaluate uncertainties on tt modelling, were generated using the Powheg- 
Box and MC@NEO v4.1 [54-56] generators interfaced with Herwig v6.5 [57, 58] with Jimmy [59] for the 
modelling of the underlying event. Eor these samples the CTIO PDE set and the ATEAS-AUET2 [60] tune 
were used. In all cases, the top quark mass used was 172.5 GeV. The cross-sections for these samples were 
normalised to the calculation from Top-p-i- v2.0 [61] at next-to-next-to-leading-order (NNEO) accuracy in 
the strong coupling constant as. including resummation of next-to-next-to-leading-logarithmic (NNEE) 
soft gluon terms [62-67]. The top quark kinematics in all SM tt samples were corrected to account for 
electroweak higher-order effects [68]. This correction is applied after generating the samples, by applying 
a weight that depends on the flavour of the initial-partons, the centre-of-mass energy of the initial partons, 
and the decay angle of the tops in the centre-of-mass frame of the initial partons. 


^ The 4% loss here is relative to the full 20.3 fb 'data set. 
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Production of W bosons in association with jets (fP+jets) is also a significant background process. 
Samples of IP+jets events were generated using the Alpgen v2. 13 generator [69] interfaced with Py- 
THiA v6.426, including up to five exfra parfons in fhe mafrix elemenf. Configurations wifh additional 
heavy quarks (a single c-quark, a cc pair or a bb pair in fhe hard process) were included, wifh fhe masses 
of heavy quarks faken info accounf. The CTEQ6L1 PDF sef and fhe Perugia 2011C fune were used. 
The samples were normalised using dafa as described in section 7. Addifional samples were generafed 
wifh differenf choices of Alpgen mafrix elemenf/parfon shower mafching paramefers in order fo esfimafe 
modelling uncerfainfies in fhe producfion of IT+Jefs evenfs. 

Production of single lop quarks can yield evenfs lhal safisfy fhe analysis even! selection. The Powheg- 
Box generafor inferfaced wifh Pythia v6.246 was used lo esfimafe fhe s - and ITf-channels [70-72] wifh 
the same configuration as for the tt samples. The t-channel was also generated with Powheg-Box but 
in a four-flavour scheme, hence fhe CTIO NLO four-flavour PDF sef was used. Overlap befween fhe 
Wt sample and ti samples was handled using fhe diagram removal scheme [73]. These samples were 
normalised fo approximale NNFO cross-sections [74-76]. 

Olher minor background processes producing prompl isolafed leptons include heavy diboson producfion, 
producfion of Z bosons in association wifh jels (Z-f-jels) and producfion of heavy gauge bosons in asso- 
ciafion wifh ti (ttV). Producfion of Z-i-jels was modelled using Alpgen inferfaced wifh Pythia v6.426, in 
fhe same contiguralion used for fhe IT-i-jels samples described above. The samples were normalised fo 
the inclusive Z boson production cross-section calculated at NNFO in QCD using Fewz [77]. Diboson 
production was modelled using the Sherpa [78-81] generator, with up to three extra partons in the matrix 
element and taking into account the mass of the b- and c-quarks. The diboson samples were normalised 
to calculations at NFO in QCD performed using Mcfm [82] v5.8. The tW production was modelled using 
MadGraphS interfaced with Pythia v6.426 and normalised to NFO cross-section predictions [83]. 


5 Event selection 

Events consistent with tt decaying to a single charged lepton together with hadronic jets and missing 
transverse momentum are selected. Electron candidates are required to have a transverse energy Fj > 
25 GeV and |?7ciusterl < 2.47, where lyduster is the pseudorapidity of the cluster of energy deposited in the 
electromagnetic calorimeter, computed with respect to the centre of ATFAS detector and matched to the 
candidate [84]. Electron candidates in the transition region 1.37 < |?7ciusterl < 1-52 between calorimeter 
barrel and endcap are excluded. 

In order to reduce backgrounds from non-prompt sources and hadronic showers with a high electromag¬ 
netic energy fraction in the calorimeter, electron isolation is imposed using a mini-isolation variable M/io 
defined as fhe ^ for all fracks (excepf fhe mafched lepfon frack) wifh pj > \ GeV satisfying 

tracks 

qualify selection criferia, and wifhin a cone of size A/? < 10 GeV/F t [16,85] cenfred on fhe barycenfre of 
fhe clusfer. Elecfrons are defined fo be isolafed if M/io/Ft < 0.05. The Ft here is fhe fransverse energy 
of fhe reconsfrucfed electron. The isolation variable M1\q is particularly useful in the case of boosted top 
quark decays since the pT-dependent cone size reflects the pT dependence of the separation between ob¬ 
jects with a common boosted parent. The matching of the electron to the collision vertex [86] is imposed 
by requiring that the longitudinal impact parameter relative to it be less than 2 mm. 

Muon candidates are required to have pT > 25 GeV and \ r ]\ < 2.5. The matching of the muon to the 
collision vertex is imposed by the requirements that the longitudinal impact parameter relative to the 
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collision vertex be less than 2 mm and that the transverse impact parameter relative to the collision vertex 
divided by its uncertainty, |r/o/o‘rfol’ ’^^an 3.0. Muons are also required to satisfy the same MI\q 

requirement as electrons, with the cone centred on the inner-detector track associated with the muon. 

Jets are reconstructed using the anti-k, algorithm [87] applied to clusters of calorimeter cells that are topo¬ 
logically connected and calibrated to the hadronic energy scale [88] using a local calibration scheme [89]. 
Small-radius jets (radius parameter R - 0.4) as well as large-radius jets {R - 1.0) are used. The ener¬ 
gies of all jets and the masses of the large-radius jets have been calibrated to their values at particle 
level [90-92]. Small-radius jets are required to satisfy pj > 25 GeV and \r]\ < 2.5, while large-radius 
jets are required to satisfy pj > 300 GeV and \r]\ < 2.0. Low-px central small-radius jets (pj < 50 GeV, 
It/I < 2.4) are required to have a jet vertex fraction [93] greater than 0.5. The jet vertex fraction is defined 
as the total transverse momentum (using a scalar sum) of tracks in the jet that are associated with the 
primary vertex divided by the scalar sum of the transverse momentum of all tracks in the jet. This vari¬ 
able suppresses jets arising from pile-up effects. Large-radius jets have jet trimming [17] applied. In 
trimming, subjets are formed by applying a jet algorithm with smaller radius parameter, /?sub> and then 
soft subjets with less than a certain fraction, /cut, of the original jet pj are removed. The properties of the 
trimmed jet are then calculated using the surviving subjets. This procedure mitigates the effect of pile- 
up [94]. The trimming parameters used in this search are /cut - 0.05 and /?sub = 0.3, and the inclusive 
kt [95] algorithm is used to form the subjets. 

Only small-radius jets are considered for fr-jet identification (fr-tagging). The fr-tagging algorithm uses a 
multivariate approach with inputs taken from the results of separate impact parameter, secondary vertex 
and decay topology algorithms [96]. The operating point of the algorithm is chosen such that the b- 
tagging efficiency for simulafed tt evenfs is 70%. In MC simulafion, facfors are applied fo correcf for fhe 
differences befween fhe fr-tagging efficiency in simulafed evenfs and fhaf measured in dafa. The facfors 
are adapfed fo be appropriafe for Zr-jefs from high-px top quarks, for which fhe fr-tagging efficiencies are 
lower. 

The is calculafed from fhe vector sum of fhe fransverse energy of topological clusfers in fhe calori¬ 
meter [97]. The clusfers associafed wifh fhe reconsfrucfed elecfrons and small-radius jefs are replaced by 
fhe calibrafed energies of fhese objecfs. Muon fransverse momenfa defermined from fhe ID and fhe muon 
specfromefer are also included in fhe calculation. 

Overlap in idenlificalion of fhe relevanf physics objecfs is possible and a procedure is implemenfed fo 
remove duplication. Elecfrons and small-radius jefs are considered for overlap removal if fhe cluster 
associafed wifh fhe elecfron is wifhin AT? = 0.4 of fhe nearesf jef. In such cases, fhe jefs have fheir four- 
momenfum and jef vertex fraclion recalculated afler subfracfing fhe election four-momenlum and fhen 
are removed if fhe recalculated values do nol safisfy fhe original jef selecfion criteria. If fhe disfance 
AT? befween fhe elecfron and fhe recalculated jef is < 0.2, fhe elecfron candidate is likely to be from fhe 
hadronic jef. Therefore fhe election is removed from fhe elecfron candidafe lisf and ifs four-momenfum 
is added fo fhaf of fhe recalculafed jef. Muons are removed from fhe muon candidafe lisf if fhe disfance 
AT? befween fhe muon and small-radius jef is less fhan 0.04 -i- 10 GeV/px- This criterion exploifs fhe 
anfi-correlalion befween fhe muon px and ifs angular disfance from fhe Tr-quark, in an approach similar 
fo fhe isolafion variable. The paramefers are funed based on signal MC simulafions in order fo provide a 
consfanf high efficiency as a funcfion of resonance mass. 

A sef of common preselection criteria is used for evenfs fo be considered for fhe boosfed and resolved 
fopologies. Evenfs are required fo have exacfly one lepton (election or muon) plus multiple jefs. Evenfs 
recorded by fhe lepton friggers are required fo have lepfon-tiigger objecfs fhaf mafch fhe selecfed lepton. 
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Additionally, events must have > 20 GeV and + wit > 60 GeV, where mj is the transverse 
mass calculated as nix = - cos(p{y) where pj is the transverse momentum of the lepton and 

(pev is the angle between the pj and vectors. 

Events are next checked against the boosted-topology selection. The selected lepton is required to have at 
least one small-radius jet within a distance of AR({, j) = 1.5 and, of these, the jet with highest px is termed 
y'sei- Boosted-topology events must have at least one large-radius jet with pj > 300 GeV (380 GeV for 
the muon-plus-jets events selected by the large-radius-jet trigger), \q\ < 2.0, mass mjet > 100 GeV, first k, 
splitting scale [95] ^fdn > 40 GeV, A/? >1.5 between the large-radius jet and jsei, and A0 > 2.3 between 
the large-radius jet and lepton. The jet mass is calculated using the four-momenta of its constituent 
clusters, which are taken as massless. If multiple large-radius jets satisfy these criteria, the highest-px jet 
is chosen as the hadronically decaying top quark candidate. Finally, at least one of the small-radius jets in 
the event must be f^-tagged and matched to either of the top quark candidates, as described in section 6. 

Events that do not satisfy the boosted-topology selection are then tested against the criteria for the 
resolved-topology selection. These events are required to have at least four small-radius jets, with at 
least one of them fj-tagged. A algorithm is used to reconstruct the tt system, as described in section 6. 
The lowestvalue is required to satisfy logjQ Of") < 0.9. 

The selection efficiencies for MC simulafed signal evenfs are given in figure 1, for differenl models of 
inferesf. For reference, fhe branching rafio for tt fo elecfron- or muon-plus-jefs is abouf 17% for each 
lepfon flavour faking info accounf lepfonic fau decays [98]. There are efficiency losses from bofh fhe 
large-radius jef requiremenfs and fj-fagging requiremenf for fhe boosfed-fopology selection. For resonance 
masses above 1.5 TeV, fhe efficiencies of fhe resolved selections are relatively insignificanf due fo fhe x^ 
requiremenfs and fhe vefo of boosfed selections. If can also be seen fhaf efficiency times accepfance 
is smaller for isolafed elecfrons fhan isolafed muons above fhe same resonance mass poinf, due fo fhe 
inefficiency of elecfron idenfificafion and overlap removal in fhe boosfed environment. 


6 Event reconstruction 

Signal ti resonances should appear in the spectrum as an excess of events over the SM expectation 
clustered around the resonance mass. Events are reconstructed assuming the final state originated from a 
tt decay. To calculate the neutrino four-momentum must be determined. The neutrino transverse 
momentum is taken to be the vector. The longitudinal component of the neutrino momentum, p^, 
is calculated by constraining the lepton plus missing momentum system to have the W boson mass and 
solving the resulting quadratic equation in the neutrino’s longitudinal momentum p^ [99, 100]. If no real 
solution exists, the F™®® vector is varied by the minimal amount required to produce exactly one real 
solution. If two real solutions are found, the one with the smallest Ip^l is used for the boosted-topology 
reconstruction, while the choice is made by the;y^ algorithm described below for the resolved topology. 

For the boosted topology, is computed from the four-momenta of the neutrino, lepton, the previously 

selected small-radius jet, and the highest-px large-radius jet. In this case the assignment of jets to the 

semileptonically decaying top quark and hadronically decaying top quark is unambiguous. 

In calculating for the resolved topology, a algorithm is employed to find the best assignment 
of jets to the semileptonically and hadronically decaying top quarks. Using the four-momenta of the 







Figure 1: Selection efficiency times acceptance times branching ratio as a function of the top-antitop quark invariant 
mass m,j at Monte Carlo generator level for the different signals in the models considered; (a) Z', (b) ^kk, (c) 
Gkk, (d) scalar resonance. The dashed lines show the boosted-topology selection and the unbroken lines show the 
combined selection. 


neutrino, lepton and all small-radius jets in the event, is defined using the expected top quark and W 
boson masses: 


nijj — 
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2 

bbl jiv 

(T w 


cr 




(PTJjb - PTjev) - (PTA ~ PT.te) ^ 
O' diffpx 


( 1 ) 


The first term is a constraint using the mass of the hadronically decaying W boson. The second term is 
a constraint using the mass difference between the hadronically decaying top quark and the hadronically 
decaying W boson. Since the mass of the hadronically decaying W boson, nijj, and the mass of the 
hadronically decaying top quark, nijj},, are highly correlated, the mass of the hadronically decaying W 
boson is subtracted from the second term so as to decouple it from the first term. The third term is a 
constraint using the mass of the semileptonically decaying top quark. The last term arises as a constraint 
on the expected transverse momentum balance between the two decaying top quarks. In the definition 
above, th and tf refer to the hadronically and semileptonically decaying top quarks. The values of the 
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central-value parameters m-^, ^tt and pj.fh ~ PT,te-> and the values of the width parameters 

cr^y, CTtt and (Ti\spj are found through Gaussian fits to the distributions of relevant reconstructed 

variables, using reconstructible^ MC events with Z' masses from 0.5 to 2.0 TeV. All possible neutrino 
solutions and jet permutations with at least one f^-quark candidate satisfying the f7-tagging requirement 
are tested, and the one with the lowest value is adopted. About 80% of these reconstructible events 
have the correct assignment of reconstructed objects to the hadronically and semileptonically decaying 
top quarks. 

The resulting distributions for several signal masses from the resolved- and boosted-topology recon¬ 
struction are shown in figures 2 and 3. For fhese figures, all evenfs safisfying fhe resolved- or boosfed- 
fopology selecfion criteria are used. The low-mass fails arise from fwo effecfs: firslly, exfra radiafion from 
the tt system that is not included in the reconstruction can shift the reconstructed mass to lower values; 
secondly, before reconstruction the Breit-Wigner signal shape in has a tail at lower values due to the 
steep fall in parton luminosity with increasing partonic centre-of-mass energy. The former is particularly 
true for high-mass resonances, while the latter has a larger effect on broad resonances. The experimental 
resolution for the invariant mass of the tt system^ is 8% for the resolved-topology selection at a resonance 
mass of 400 GeV improving to 6% for 1 TeV. It is 6% in the boosted-topology selection, independent of 
resonance mass. 

With hadronically and semileptonically decaying top quarks identified for bofh fhe boosfed- and resolved- 
fopology selections, fhree cafegories of f^-tagged evenfs are defined: fhose in which bofh decaying fop 
quark candidates have a mafching b-jet, fhose in which only fhe hadronically decaying lop quark candidale 
has a mafching f^-jel and fhose in which only fhe semilepfonically decaying lop quark candidate has a 
mafching fj-jel. In fhe boosled-lopology selection, fhe malch is defined as eilher fhe selected jel, jseh or a 
small-radius jel wilhin AR = 1.0 of fhe large-radius jel as being a b-tagged jel. In fhe resolved-fopology 
selection, fhe mafching is determined by fhe x^ algorifhm. 


7 Background contributions estimated from data 


Dafa are used lo eslimafe fhe magnifudes and uncerlainlies of fwo imporlanf background conlribulions: 
W-rjels and multi-jel production. 


7.1 W-i-jets scale factors 

For fhe IT-tjels background, dafa are used lo derive scale faclors lhal are applied lo correcf fhe normalisa¬ 
tion and flavour fracfions given by Alpgen MC simulafions of fhis background. 

For bofh fhe resolved- and boosled-lopology evenl selecfion crileria, fhe normalisalion scale factors are 
delermined by comparing fhe measured W boson charge asymmelry in dafa [101,102] wilh lhaf predicled 
by Alpgen MC simulafion. For fhe resolved fopology, all selecfion crileria are applied fo fhe dafa excepf 

^ Reconstructible events are those where there is a reconstructed object within AR = 0.4 of each visible parton from top decays, 
and the Acp between the neutrino and the is smaller than 1. Among the events satisfying the final selection critera, 70% 
(55%) are reconstructible for Z' mass of 750 GeV (400 GeV). 

^reco,true—matched _ 

* The experimental resolution is extracted from a Gaussian fit of the relative difference — - - from reconstructible 

events; reconstructed mT“ computed with the correct combination of jets identified with the parton-level 

information and m,f is the true mass of the tt system in the MC simulation. 
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Figure 2: Reconstructed top quark pair invariant mass, for the different signal models for events satisfying 
the resolved-topology selection and using the resolved-topology reconstruction: (a) Z', (b) 0 kk. (c) Gkk. (d) scalar 
resonance, before and after the veto on the boosted selection. For each signal, the two histograms are normalised to 
the same arbitrary luminosity. 


for the Zt-tagging requirement. For the boosted topology, in order to decrease the statistical uncertainty on 
the scale factors, a relaxed set of selection criteria that does not include the Z>-tagging, A0(jet, 1) > 2.3, jet 
mass and requirements is used. Any bias induced by relaxing the criteria for the boosted selection is 
found to be negligible compared to the statistical uncertainty in the scale factor determination. The total 
number of IT-i-jets events in data, Nw+ + Nw-, is given by: 

Nw^ +Nw-=i (Dcon-+ - £>con-), (2) 

Vue - 1 / 

where rue is the ratio given by MC simulation of the number of IT-i-jets events with a positively charged 
lepton to that with a negatively charged lepton and Dcon-+(-) is the number of observed events with a 
positively (negatively) charged lepton. Contributions to Dcorr+(-) from charge-asymmetric processes such 
as single top, ITZ and tt+W boson production are estimated from MC simulation and are subtracted. 
Contributions from charge-symmetric processes such as tt production cancel in the difference on the 
right-hand side of eq. (2). A scale factor, Ca, applied to the MC simulated samples of IT -i- jets events 
is then calculated as the ratio of Nw+ + Nw- evaluated from data to that predicted from MC simulation. 
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Figure 3; Reconstructed tt invariant mass for the different signal models for events satisfying the boosted-topology 
selection and using the boosted-topology reconstruction; (a) Z', (b) 0 kk. (c) Gkk. (d) scalar resonance. 


The value and statistical uncertainty obtained for Ca in the electron (muon) channel are 1.026 + 0.011 
(0.978 ± 0.010) with the resolved selection, and 0.89 ± 0.06 (0.81 + 0.05) with the boosted selection. 

Scale factors for the relative fraction of heavy-flavour contributions from W + bb, W + cc, W + c 
also determined from data [102-104]. In determining these scale factors, events are required to satisfy 
all selection criteria common to the boosted- and resolved- topology selections. Exactly two small-radius 
jets are required without any tagging requirement. 

The flavour fractions for W + bb, W + cc,W + c and IT-i-light-quark flavours are first determined from MC 
simulation for this sample. The ratio of the W + bb to W + cc contribution is taken from MC simulation 
and fixed at that value. A system of three equations is used to fit to the two-jet data sample with at least 
one (7-tagged jet in order to determine correction factors for each of the flavour fractions determined from 
MC simulation: 
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where Dvr± is the expected number of VK+jets events with a positively charged or negatively charged 
lepton in the data. The flavour fraction is /fiavour^ and the correction factor for a given flavour component 
is ^flavour- The different flavour labels are bb, cc, c, and light corresponding to IT + bb, W + cc,W + c, and 
IT+light-jets respectively. The numbers of positively charged and negatively charged leptons in the data 
are found by subtracting all non-lT+jets contributions, which are determined from MC simulations as 
35% (15%) of the selected events for the electron (muon) channel. An iterative process is used to find the 
^flavour factors, which are then used to correct the corresponding flavour fractions /flavour that are applied 
during the Ca factor calculation. In this iterative process, only the A'flavour and Ca factors are allowed to 
vary. The .fifflavour factors are initially set to unity, thus altering the Ca factor calculation. New correction 
factors ^'flavour are calculated by inverting eq. (3), and then the process is repeated ten times, with each 
repetition using the correction factors determined from the previous one. It was checked that using more 
than ten iterations produces only negligible changes in the extracted correction factors. 

The correction factors found for the two-jet sample are extrapolated to higher jet multiplicities by keeping 
the same ratios between them while conserving the normalisation in each jet multiplicity bin. The .fifflavour 
factors thus obtained are different from unity. For events containing electrons (muons), the extracted 
values and statistical uncertainties of and K^c are 1.36 ± 0.07 (1.51 ± 0.08), is 0.71 ± 0.03 (0.66 ± 
0.03), and T^nght is 0.934 ± 0.005 (0.873 ± 0.004). 

This method reduces the systematic uncertainties (from the jet energy scale, Z?-tagging and other uncer¬ 
tainties) compared to using VT-i-jets MC simulation alone. Systematic uncertainties in the IT-i-jets norm¬ 
alisation and flavour-fraction corrections are determined by rederiving these scale factors when a given 
systematic effect is applied. The new scale factors are then used in producing the mass spectrum for 
that particular systematic uncertainty. 


7.2 Multijet estimate 

The multi-jet background in events satisfying the resolved or boosted selection criteria consists of events 
with a jet that is misreconstructed as a lepton or with a non-prompt lepton that satisfies the identification 
criteria. The normalisation, shape, and statistical and systematic uncertainties associated with the 
multi-jet background are estimated from data using a matrix method [103,105]. 

The matrix method utilises efficiencies for leptons produced by prompt and non-prompt sources. The 
efficiency / is defined as the probability that a non-prompt lepton from multi-jet production that satisfies 
the loose identification criteria [105] also satisfies the tight identification criteria. It is derived from data 
in control regions dominated by multi-jet events, with prompt-lepton contributions subtracted based on 
MC simulations. The efficiency e is defined as the probability that a lepton from prompt sources (W or 
Z bosons) that satisfies the loose identification criteria also satisfies the tight identification criteria. It 
is determined using SM MC samples with a mixture of processes similar to that in the signal region, 
corrected using data versus MC correction factors derived from Z ^ events. 

The number of multi-jet background events satisfying the resolved or boosted selection criteria is estim¬ 
ated using data events that satisfy all selection criteria, except that the loose lepton identification criteria 
are used. This sample contains prompt as well as non-prompt leptons. 

The number of events with leptons satisfying the loose identification criteria, Al is defined as 

— ^prompt + ^multi-jet (4) 
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where A^piompt is the number of events with prompt leptons satisfying the loose identification criteria and 
iVmuiti-jet is the number of events satisfying the loose identification criteria with leptons from other sources. 
The number of events satisfying the tight identification criteria, is then 

f X A^prompt X A^multi-jet • (5) 

Solving these two equations for A^prompt and A^muiti-jet gives the multi-jet contribution from events satisfying 
all the selection criteria. 

Good shape modelling of the distributions is achieved by parameterising the efficiencies as functions 
of relevant kinematic variables, and validated in the multi-jet control regions. Systematic uncertainties 
are evaluated using several different definitions of multi-jet control regions that result in slightly different 
/ estimations. Systematic uncertainties associated with object reconstruction and MC simulation are also 
considered, resulting in a total normalisation uncertainty of 20%. 


8 Systematic uncertainties 


The systematic uncertainties can be broadly divided into two categories: uncertainties that affect recon¬ 
structed objects (such as jets) and uncertainties that affect the modelling of certain background or signal 
processes. Some of the uncertainties affect both the shape and the normalisation of the spectra, 
while others affect the normalisation only. In table 1 , an overview of the effects of the dominant system¬ 
atic uncertainties on the background and signal yields is given. Only the impact on the overall normal¬ 
isation is shown in the table, but some of the systematic uncertainties have a significant dependence on 
the reconstructed tt mass, which is fully taken into account in the analysis for all of these uncertainties. 
The systematic uncertainties with the strongest dependence are those from the jet energy scale, parton 
distribution functions and A-jet identification. 

The dominant uncertainty on the normalisation of the total background estimate is the NNLO-i-NNLL 
tt cross-section uncertainty of 6.5%. This uncertainty includes renormalisation and factorisation scale 
uncertainties, combined PDF and strong-coupling uncertainties evaluated following the PDF4LHC [106] 
recommendations, and uncertainties associated with the value of the top quark mass. The combined PDF 
and strong-coupling uncertainties are extracted for each of the three PDF sets: MSTW2008 68% confid¬ 
ence level (CL) NNLO [43,107], CTIO NNLO [52,108] and NNPDF2.3 NNLO [109]; the total uncer¬ 
tainty associated with the PDF and strong-coupling uncertainties is one half of the size of the envelope of 
the three resultant error bands, with the central prediction being the midpoint of the envelope. Variations 
from changing the top quark mass by ±1.0 GeV are added in quadrature to the scale uncertainties, and 
combined PDF and strong-coupling uncertainties. 

The VT±jets normalisation, determined using data, is separately evaluated for each experimental source of 
systematic uncertainty. An additional systematic uncertainty on the prediction is evaluated by using the 
simulated samples generated with varied Alpgen matching parameters. The normalisation uncertainty on 
the multi-jet background is 20%, as described in section 7. 

The single-top-quark background normalisation uncertainty is 7.7% [74-76]. The normalisation uncer¬ 
tainty on the Z±jets sample is 48%, estimated using Berends-Giele scaling [110]. The diboson normalisa¬ 
tion uncertainty is 34%, which is a combination of the NLO PDF and scale uncertainties, and additional 
uncertainties from the requirements on the jet multiplicity. 
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Table 1: Average impact of the dominant systematic uncertainties on the total background yield and on the estim¬ 
ated yield for a Z' sample with m = 1.75 TeV. The electron and muon channel spectra are added. Any shift in 
yield is given in percent of the nominal value. Certain systematic uncertainties are not applicable to the Z' samples, 
which is indicated with a bar (-) in the table. 



Resolved selection 

Boosted selection 


yield impact [%] 

yield impact [%] 

Systematic Uncertainties 

total bkg. 

Z' 

total bkg. 

Z' 

Luminosity 

2.5 

2.8 

2.6 

2.8 

PDF 

2.4 

3.6 

4.7 

2.3 

ISR/FSR 

3.7 

- 

1.2 

- 

Parton shower and fragmentation 

4.8 

- 

1.5 

- 

tt normalisation 

5.3 

- 

5.5 

- 

tt EW virtual correction 

0.2 

- 

0.5 

- 

tt generator 

0.3 

- 

2.6 

- 

tt top quark mass 

0.6 

- 

1.4 

- 

fk-l-jets generator 

0.3 

- 

0.1 

- 

Multi-jet normalisation, e-l-jets 

0.5 

- 

0.2 

- 

Multi-jet normalisation, yu+jets 

0.1 

- 

<0.1 

- 

JES-l-JMS, large-radius jets 

0.1 

2.1 

9.7 

2.8 

JER+JMR, large-radius jets 

<0.1 

0.3 

1.0 

0.2 

JES, small-radius jets 

5.6 

2.6 

0.4 

1.4 

JER, small-radius jets 

1.8 

1.4 

<0.1 

0.2 

Jet vertex fraction 

0.8 

0.8 

0.2 

<0.1 

h-tagging J>-jet efficiency 

1.1 

2.0 

2.9 

17.1 

h-tagging c-jet efficiency 

0.1 

0.7 

0.1 

2.1 

f>-tagging light-jet efficiency 

<0.1 

< 0.1 

0.5 

0.2 

Electron efficiency 

0.3 

0.6 

0.6 

1.3 

Muon efficiency 

0.9 

1.0 

1.0 

1.1 

MC statistical uncertainty 

0.4 

6.0 

1.3 

1.8 

All systematic uncertainties 

10.8 

8.8 

13.4 

18.0 


The uncertainty on the integrated luminosity is 2.8%. It is derived, following the same methodology as 
that detailed in Ref. [Ill], from a preliminary calibration of the luminosity scale derived from beam- 
separation scans performed in November 2012. This uncertainty is applied to all signal and background 
samples except multi-jet and VT-i-jets, which are estimated from data. 

The effect of the PDF uncertainty on all MC samples is estimated by taking the envelope of the NNPDF2.3, 
MSTW2008NLO and CTIO PDF set uncertainties at 68% CL^ following the PDF4LHC recommendation 
and normalising to the nominal cross-section. The PDF uncertainty on the tt mass spectrum has a much 
larger effect on the boosted sample than on the resolved sample. The effect on the total background yield 
is 2.1% (4.2%) after the resolved (boosted) selection. The size of the uncertainty grows with reconstructed 
mass, attaining values of 50% above 2 TeV in the boosted selection. 

One of the dominant uncertainties affecting reconstructed objects is the jet energy scale (JES) uncertainty, 
especially for large-radius jets [91,92]. Uncertainties on the jet mass scale (JMS) and the kt splitting 
scales [92] are also important for this analysis. These uncertainties have an impact of 10% on the overall 
background yield in the boosted selection. The impact on the background estimates falls with increas¬ 
ing varying from 22% at lowest masses to about 7% above 1.5 TeV. The impact is smaller for the 

^ The CTIO PDF uncertainties are scaled down by a factor 1.645 to reach an approximate 68% CL. 
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resolved selection, since the large-radius jets are only used indirectly there via the vetoing of events that 
satisfy the boosted selection. For large-radius jets, uncertainties on the jet energy resolution (JER) and 
jet mass resolution (JMR) are also considered. These are less significant than the JES. Eor small-radius 
jets, the uncertainties on the JES, the jet reconstruction efficiency and the jet energy resolution are con¬ 
sidered [88,90]. The small-radius JES uncertainty is one of the most significant systematic uncertainties 
in the resolved selection, affecting the overall expected yield by 6%. The effect of uncertainties associated 
with the jet vertex fraction is also considered. The f^-tagging uncertainty is modelled through simultan¬ 
eous variations of the uncertainties on the efficiency and rejection [112, 113]. An additional fj-tagging 
uncertainty^ is applied for high-momentum jets {pi > 300 GeV) to account for uncertainties on the mod¬ 
elling of the track reconstruction in high-pj environments, which is one of the dominant uncertainties for 
high-mass signals. 

Eor the leptons, the uncertainties on the isolation efficiency, the single-lepton trigger and the reconstruc¬ 
tion efficiency are estimated using Z ^ ee and Z ^ pp events. In addition, high-jet-multiplicity Z ^ 
events are studied, from which extra uncertainties on the isolation efficiency are assigned to account for 
the difference befween Z and tt evenfs. Uncerfainfies on fhe reconsfrucfion, as well as on fhe energy 
scale and resolufion of fhe leptons, are also considered, and generally have a smaller impacf on fhe yield 
and fhe search sensifivify fhan fhe uncerfainfies menfioned above. 

The uncerfainfy on fhe tt background due fo uncerfainfies on fhe modelling of QCD initial- and final- 
sfafe radiation (ISR/ESR) is estimated using AcerMC v3.8 [114] plus Pythia v6.426 MC samples by 
varying fhe Pythia ISR and ESR parameters wifhin ranges allowed by a previous ATEAS measuremenf 
of tt production wifh a vefo on addifional cenfral jef acfivify [115]. The QED ISR/ESR uncerfainly is 
negligible af fhis level. The dependency of fhe shape on fhe choice of NEO generafor is accounfed 
for by using fhe difference befween samples generafed wifh MC@NEO and Powheg-Box-i-Herwig as 
a systematic uncerfainfy. The parfon showering and fragmenfafion uncerfainfy on fhe tt background is 
esfimafed by comparing fhe resulfs from a sample generated wifh Powheg-Box [48], interfaced wifh 
Pythia or Herwig. The uncerfainly on fhe dislribulion arising from fhe uncertainly in fhe fop quark 
mass is evaluated by comparing fhe speclrum using fhe nominal sample fo Ihose generafed wifh lop 
quark masses of 170 and 175 GeV, and multiplying fhe difference by 0.4 (fo approximale a one slandard 
devialion uncertainly, corresponding fo ±1.0 GeV). The uncertainly on fhe eleclroweak correcfions fo 
top quark pair production is modelled by changing fhe difference of each correction factor from unify by 
±10%. 

Eor fhe lT±jefs background, fhe uncertainly on fhe disfribufion is esfimafed by reweighling fhe 
evenfs fo fhe kinemalics of MC samples generafed wifh a differenl malching scale or functional form of 
fhe facforisafion scale [69]. 


9 Comparison of data with expected background contributions 


Afler all evenl selection criteria are applied, 223,330 evenfs remain for fhe resolved lopology and 8,206 
evenfs are selecfed for fhe boosfed lopology. The evenl yields from dala and from expecled background 
processes are listed in fable 2 fogefher wifh fhe associaled syslemalic uncerfainfies. 


^ The additional tagging uncertainty is an extrapolation of the uncertainty from regions of lower pj. Depending on the px of 
the jet, it is approximately 12% to 33% for h-jets and 17% to 30% for c-jets, added in quadrature to the uncertainty on the jet 
h-tagging efficiency correction factor for the 200-300 GeV region. 
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Table 2: Data and expected background event yields after the resolved and boosted selections. The sum in quadrat¬ 
ure of all systematic uncertainties on the expected background yields is also given. 



Resolved-fopology selecfion 


Type 

e-i-jels 

fi+jets 

Sum 

tt 

93,000 ± 11,000 

91,000 ± 11,000 

184,000 ± 22,000 

Single fop 

3,800 ± 500 

3,800 ± 500 

7,600 ± 1,000 

ttV 

274 ± 40 

267 ± 40 

541 ± 80 

Mulfi-jef e 

5,300 ± 1,100 

- 

5,300 ± 1,100 

Mulfi-jef ju 

- 

1,050 ± 240 

1,050 ±240 

IF-i-jefs 

6,600 ± 800 

7,100 ± 800 

13,700 ± 1,500 

Z-tjefs 

1,400 ± 750 

650 ± 340 

2,000 ± 1,080 

Dibosons 

320 ± 120 

310 ± 120 

620 ± 240 

Tofal 

110,000 ± 12,000 

105,000 ± 12,000 

215,000 ± 24,000 

Dafa 

114,377 

108,953 

223,330 


Boosfed-fopology selection 


Type 

e-i-jels 

//-i-jels 

Sum 

tt 

4,100 ± 600 

4,000 ± 600 

8,100 ± 1,200 

Single fop 

138 ± 20 

154 ± 20 

290 ± 40 

ttV 

37 ±6 

38 ±7 

75 ± 13 

Mulfi-jef e 

91 ± 18 

- 

91 ± 18 

Mulfi-jef n 

- 

8.6 ± 1.6 

8.6 ± 1.6 

VF-i-jefs 

260 ± 50 

290 ± 50 

550 ± 100 

Z-i-jefs 

31 ± 16 

17 ± 9 

48 ±25 

Dibosons 

21 ± 8 

20 ± 8 

41 ± 16 

Tofal 

4,700 ± 600 

4,500 ± 600 

9,200 ± 1,200 

Dafa 

4,148 

4,058 

8,206 


Good agreement is observed between the data and the total expected background. Figure 4 shows the 
reconstructed mass of the semileptonically and hadronically decaying top quark candidates, and the mass 
of the hadronically decaying W boson candidate for the resolved selection. The equivalent distributions 
for the boosted selection are the mass of the semileptonically decaying top quark candidate and the mass 
of the large-radius jet; both are shown in figure 5. Figure 6 shows the distribution of the transverse 
momentum and first k; splitting scale of the selected large-radius jets. 

The tt invariant mass spectra for the resolved and the boosted selections in the electron and muon chan¬ 
nels, separated by f7-tagging category, are shown in figures 7 and 8. Figure 9 shows fhe tt invarianf mass 
specfrum for all channels added fogefher. The dafa generally agree wifh fhe expecfed background wifh 
slighf shape differences seen especially in fhe high-mass regions; fhese deviafions are consisfenf wifh fhe 
uncerfainfies. 
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(a) e+jets channel. 


(b) yu+jets channel. 
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(c) e+jets channel. 


(d) ^+jets channel. 




(e) e+jets channel. 


(f) yu+jets channel. 


Figure 4: Reconstructed mass of the hadronically decaying top quark candidate, semileptonically decaying top 
quark candidate, and hadronically decaying Vk-boson candidate after the resolved-topology selection in the electron 
and muon channels. The SM background components are shown as stacked histograms. The shaded areas indicate 
the total systematic uncertainties. 
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(a) e+jets channel. 


(b) yu+jets channel. 



100 120 140 160 180 200 220 240 260 280 300 
Semileptonic top mass [GeVJ 



Semileptonic top mass [GeVJ 


(c) e+jets channel. 


(d) ^+jets channel. 


Figure 5: The invariant mass of the large-radius jets and the invariant mass of the semileptonically decaying top 
quark candidate, after the boosted selection. The SM background components are shown as stacked histograms. 
The shaded areas indicate the total systematic uncertainties. 
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(a) e-l-jets channel. 


(b) yu-l-jets channel. 
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(c) e-l-jets channel. 


(d) yu-l-jets channel. 


Figure 6: The transverse momentum, pj , and first k, splitting scale, yfdn, of the large-radius jet after the boosted 
selection. The SM background components are shown as stacked histograms. The shaded areas indicate the total 
systematic uncertainties. 
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(a) Electron channel, resolved selection, both 
top candidates have h-tagged jets. 



(c) Electron channel, resolved selection, only 
the hadronic top candidate has a h-tagged jet. 



(e) Electron channel, resolved selection, only 
the semileptonic top candidate has a fc-tagged jet. 



(b) Muon channel, resolved selection, both top 
candidates have fo-tagged jets. 



(d) Muon channel, resolved selection, only the 
hadronic top candidate has a i-tagged jet. 



(f) Muon channel, resolved selection, only the 
semileptonic top candidate has a Z^-tagged jet. 


Figure 7: The spectrum of the reconstructed ttbar invariant mass, for the different channels, before any nuis¬ 
ance parameter fit, after the resolved-topology selection. The SM background components are shown as stacked 
histograms. The shaded areas indicate the total systematic uncertainties. The green line shows the expected distri¬ 
bution for a hypothetical gfKK of mass 0.8 TeV, width 15.3%. 
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(a) Electron channel, boosted selection, both top 
candidates have fe-tagged jets. 



(c) Electron channel, boosted selection, only the 
hadronic top candidate has a Zj-tagged jet. 



(e) Electron channel, boosted selection, only the 
semileptonic top candidate has a b-tagged jet. 



(b) Muon channel, boosted selection, both top 
candidates have fo-tagged jets. 



(d) Muon channel, boosted selection, only the 
hadronic top candidate has a i-tagged jet. 



(f) Muon channel, boosted selection, only the 
semileptonic top candidate has a Z>-tagged jet. 


Figure 8; The spectrum for the different channels, before any nuisance parameter fit, after the boosted- 
topology selection. The SM background components are shown as stacked histograms. The shaded areas indicate 
the total systematic uncertainties. The red line shows the expected distribution for a hypothetical qvv of mass 2.0 
TeV, width 15.3%. 
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(a) Boosted selections. 


(b) Resolved selections. 



(c) All selections. 


Figure 9: The distributions, before any nuisance parameter fit, summed over (a) all 6 boosted channels, (b) all 
6 resolved channels, and (c) all 12 channels compared with data. The SM background components are shown as 
stacked histograms. The shaded areas indicate the total systematic uncertainties. The red (green) line shows the 
expected distribution for a hypothetical 0 kk of mass 2.0 (0.8) TeV, width 15.3%. 
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10 Results 


The final discriminating observables that are used to search for a massive resonance are the twelve ti 
invariant mass spectra: three b-i&g categories for two selections and two decay channels. After the re¬ 
construction of the tt mass spectra, the data and expected background distributions are compared using 
BumpHunter [116], which is a hypothesis-testing tool that searches the data for local excesses or deficits 
compared to the expected background, taking the look-elsewhere effect [117] into account over the full 
mass spectrum. The search is performed on three combinations of the spectra: the six channels of the 
resolved selections, the six channels of the boosted selections and the twelve channels. The most sig¬ 
nificant deviation of the data from the expected background spectrum is required to appear at the same 
place in each of the channels of a combination. After accounting for the systematic uncertainties, no 
significant deviation from the total expected background is found. Upper limits are set on the cross- 
section times branching ratio for each of the signal models using a profile likelihood-ratio test. The CL^ 
prescription [118] is used to derive one-sided 95% CL limits. The results are obtained using the Hist- 
Fitter [119] framework with all spectra from the 12 channels, excluding bins with few events with 
below 400 GeV in the boosted channels or above 2 TeV in the resolved channels. 

The statistical and systematic uncertainties on the expected distributions are included in this CL^ proced¬ 
ure as nuisance parameters in the likelihood fits. The nuisance parameters for the systematic uncertainties 
are constrained by a Gaussian probability density function with a width corresponding to the size of the 
uncertainty considered. Correlations between different channels and bins are taken into account. The 
product of the various probability density functions forms the likelihood function that is maximised in the 
fit by adjusting the free parameter (the signal strength) and nuisance parameters. The expected dis¬ 
tributions are compared to data in figure 10 after a fit of nuisance parameters under the background-only 
hypothesis. It can be seen that the uncertainties are smaller than in figure 9 and that the procedure is able 
to produce a good-quality fit to the data. 

The expected and observed limits using different signal models are presented in figure 11. Here the 
expected limits are obtained by taking the nominal background estimates as the expected data. For the 
ZjQ 2 benchmark, limits on the production cross-sections vary from 4.2 pb to 0.03 pb for masses from 
0.4 TeV to 3 TeV. A width 1.2% is excluded for masses less than 1.8 TeV, while masses below 

2.0 TeV are expected to be excluded. The Z'^q 2 limits are stronger for a width of 2% (3%), reaching 
2.0 TeV (2.3 TeV). For the ^kk benchmark, limits on the production cross-sections vary from 4.8 pb for 
a mass of 0.4 TeV, to 0.09 pb for a mass of 3 TeV. A ^kk of width 15.3% is excluded for masses less than 
2.2 TeV, while masses below 2.3 TeV are expected to be excluded. The cross-section limits in the Gkk 
model range from 2.5 pb for a mass of 0.4 TeV to 0.03 pb for 2.5 TeV, with no mass range excluded in 
the benchmark scenario. The cross-section limits on a narrow scalar resonance range from 3.0 pb for a 
mass of 0.4 TeV to 0.03 pb for 2.5 TeV. The cross-section limits are generally stronger for the latter two 
benchmark models than those for the spin-1 resonances due to higher acceptance. 

The width dependence of the cross-section limits was also evaluated for the ^kk models. The results are 
presented in figure 12. For a 1 TeV resonance, the limits weaken by approximately a factor of two as 
the width increases from 10% to 40%. The effect is stronger for 2 TeV and 3 TeV resonances, where the 
limits weaken by a factor of three over this width range. 

The observed upper limits on the cross-section times tt branching ratio are larger than the expected limits, 
especially for tt resonance masses greater than 1.8 TeV. This arises from the use of the profile likelihood 
method which allows the data to constrain the systematic uncertainties using the full distribution. 
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Figure 10: The distributions, after the nuisance-parameter ht under the background-only hypothesis, summed 
over (a) all 6 boosted channels, (b) all 6 resolved channels, and (c) all 12 channels compared with data. The 
SM background components are shown as stacked histograms. The shaded areas indicate the total systematic 
uncertainties. The red (green) line shows the expected distribution for a hypothetical pkk of mass 2.0 (0.8) TeV, 
width 15.3%. 
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(a) Z', resolved and boosted combination. 



(c) Gkk> resolved and boosted combination. 



(b) g'KK. resolved and boosted combination. 



.. . . 

0.5 1 1.5 2 2.5 3 


scalar resonance mass [TeV] 

(d) Scalar resonance, resolved and boosted combination. 


Figure 11; Observed and expected upper limits on the production cross-section times branching ratio to ti hnal 
states as a function of the mass of (a) Topcolour-assisted-technicolour (b) Bulk RS Kaluza-Klein gluon, 
(c) Bulk RS Kaluza-Klein graviton, (d) scalar resonance. The expected limits are derived from nominal (pre-fit) 
background estimates. The theoretical predictions for the production cross-section times branching ratio at the 
corresponding masses are also shown. 
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Figure 12: Observed and expected upper limits on the production cross-section for a Kaluza-Klein gluon times 
its branching ratio to tt, as a function of its width, for three representative mass values. The expected limits are 
derived from nominal (pre-ht) background estimates. The theoretical predictions for the production cross-section 
times branching ratio at the corresponding widths are also shown. 


thanks to the abundant data. The maximisation of the likelihood can change the central values of the 
nuisance parameters and their associated uncertainties. In the region of above 1.5 TeV, the back¬ 
ground prediction as seen in figure 9 is slightly higher than the data, which would lead to the anticipation 
that the observed limits should be slightly better than the expected ones. However, the central values of 
some nuisance parameters are significantly shifted in the fit. In the high-mj“° region one of the dominant 
uncertainties is the high-px jet f?-tagging extrapolation uncertainty, as detailed in section 8. This uncer¬ 
tainty is reduced to half of its original size by the fit, and the central value is also shifted downwards by 
approximately one pre-fit standard deviation. As this uncertainty is associated with reconstructed physics 
objects, it has correlated effects on the predictions of BSM signal and SM background. The Icr change 
in the central value of this nuisance parameter reduces the acceptance of high-mass signals in the boosted 
selection considerably: approximately 25% for a Z' with a mass of 2 TeV. As a consequence, the observed 
upper limits on the cross-section obtained from the fit to the data are larger than the expected ones fitted to 
the nominal background estimates. Constraints from the fit are also observed in the nuisance parameters 
associated with other major systematic uncertainties such as the PDF, the shape of IT-i-jets background, 
the energy and mass scales for both the large-radius and small-radius jets, f^-tagging efficiencies, and tt 
MC modelling. These constraints are understood through their impact on the fitted distributions. A 
few of the other nuisance parameters also have their central values changed slightly by the fit, but their 
impact on the signal acceptance and hence cross-section limits are much smaller than the one detailed 
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above. 


Studies were done to understand the origin of the shifted nuisanee parameter assoeiated with the high-px 
jet Z>-tagging extrapolation uneertainty, in partieular to investigate its potential origin as a mis-modelling 
of the tt baekground. Indeed an equally good fit ean be obtained by introdueing an additional ad hoc 
uneertainty in the modelling of the SM ft baekground with a similar dependeney on Sueh a mis- 

modelling would not reduee the signal aeeeptanee and results in more stringent limits on the eross-seetion 
times ft branehing ratio. Ultimately, no method was found to unambiguously aseertain the origin of this 
mis-modelling. The limits presented here are those with a Icr ehange in the fitted value of the nuisanee 
parameter assoeiated with the high-px jet f?-tagging extrapolation uneertainty, whieh eorrespond to more 
eonservative results. An upper eross-seetion limit was generated using expeeted data eonstrueted from 
a baekground-only model built from the nuisanee parameters fitted to real data. It exeludes a ^kk of 
width 15.3% for masses less than 2.1TeV, eompared to 2.2 TeV generated using nominal baekground 
estimates. 


11 Summary 


A seareh for heavy partieles deeaying to tt in the lepton-plus-jets deeay ehannel was earried out with 
the ATLAS experiment at the LHC. The seareh uses data eorresponding to an integrated luminosity of 
20.3 fb“^ of proton-proton eollisions at a eentre-of-mass energy of 8 TeV. No exeess of events beyond the 
Standard Model predietions is observed in the tt invariant mass speetra. Upper limits on the eross-seetion 
times branehing ratio are set for four different signal models: a narrow (< 3% width) Z' boson, a broad 
(15.3% width) Randall-Sundrum Kaluza-Klein gluon, a Bulk Randall-Sundrum Kaluza-Klein graviton, 
and a narrow sealar resonanee. Based on these results, the existenee of a narrow leptophobie topeolour 
Z' in the range 0.4 TeV < mz’ <1.8 TeV is exeluded at 95% CL. A broad Kaluza-Klein gluon with mass 
between 0.4 TeV and 2.2 TeV is also exeluded at 95% CL. These results probe new physies at higher mass 
than previous ATLAS searehes for the same signature, and the results are applieable to a broader variety 
of heavy resonanees. 
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